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The DEAP-1 low-background liquid argon detector has been used to measure scintillation pulse 
shapes of /3 decays and nuclear recoil events and to demonstrate the feasibility of pulse-shape discrim- 
ination down to an electron-equivalent energy of 20 keVco. The relative intensities of singlet /triplet 
states in liquid argon have been measured as a function of energy between 15 and 500 keVoo for 
both j3 and nuclear recoils. Using a triple-coincidence tag we find the fraction of /3 events that are 
mis-identified as nuclear recoils to be < 6 x 10~* between 43-86 keVcc and that the discrimination 
parameter agrees with a simple analytic model. The discrimination measurement is currently lim- 
ited by nuclear recoils induced by cosmic-ray generated neutrons, and is expected to improve by 
operating the detector underground at SNOLAB. The analytic model predicts a /3 mis-identification 
fraction of 10"^" for an electron-equivalent energy threshold of 20 keVco- This reduction allows for 
a sensitive search for spin-independent scattering of WIMPs from 1000 kg of liquid argon with a 
WIMP-nucIeon cross-section sensitivity of 10~*®cm^. 

PACS numbers: 29. 00.00,29. 25.-t,29.40.Mc 



I. INTRODUCTION 

Separation of [3-^ interactions from nuclear recoils is 
critical for many nuclear and particle astrophysics exper- 
iments, including direct searches for dark matter parti- 
cles. Liquid argon provides very sensitive pulse-shape 
discrimination based on scintillation timing and it 
is a favourable target for dark matter particle searches 
since it can be used to construct a very large target mass 
detector. In this paper we present results on the pulse- 
shape discrimination of /?-7 events from nuclear recoils 
with the DEAP-1 liquid argon detector. 

Argon has many desirable properties as a scintillator, 
among them a high light yield of approximately 40 pho- 
tons per keV [3] and ease of purification so that it can 
meet the radio-purity requirements of a rare-event search 
experiment. 
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Scintillation in argon is a result of the formation of 
dimers after exposure to ionizing radiation Q. These 
occur in singlet and triplet states. On decaying to the 
ground state, they emit light at a peak wavelength of 
128 nm, lower in energy than the lowest excited atomic 
state [1, Q . The scintillation light can thus pass through 
the argon without being absorbed. 



The two argon dimer states have vastly different life- 
timeSjEibout 6 ns for the singlet and 1.6 /xs for the triplet 
state . Moreover, the relative population of singlet and 
triplet states is determined by the linear energy trans- 
fer (LET) of the exciting radiation Q, such that fewer 
triplet excimers are produced at higher LET. With the 
large difference in lifetimes, the percentage of light signal 
in the first few tens of nanoseconds is a good estimate 
for the relative population of the singlet state, allowing 
for an effective way to discriminate between particles of 
different LETs, such as low-energy electrons and nuclei. 



II. EXPERIMENTAL APPARATUS 

The liquid argon cell consists of a stainless steel ultra- 
high vacuum chamber with two 6-inch diameter glass 
windows mounted as shown in Fig. [T] The chamber con- 
tains a volume of 5.1 L and a mass of 7 kg of liquid 
argon. All detector components have been selected to 
minimize radioactivity in the cell, in particular neutrons 
from (a,n) interactions in detector materials. A i-inch 
thick polymcthyl-mcthacrylate (PMMA) acrylic sleeve 
coated with Ti02 is inserted into the target chamber to 
act as a diffuse reflector and increase scintillation light 
collection. Scintillation photons in liquid argon have a 
mean wavelength of 128 nm. Standard vacuum deposi- 
tion techniques are used to coat the Ti02 surface of the 
acrylic sleeve and the inner surfaces of the windows with 
0.1 mg/cm^ of l,l,4,4-Tetraphenyl-l,3-butadiene (TPB) 
so as to convert the 128 nm scintillation photons into the 
visible range for detection by photomultipliers (PMTs). 
Light from the scintillation cell is transported to the 
PMTs through 8-inch long PMMA acryhc hght guides. 
To reduce potential backgrounds from Cerenkov radia- 
tion generated by cosmic-ray muons, the light guides are 
constructed of UV-absorbing (UVA) acrylic. The stain- 
less steel scintillation cell and light guides are thermally 
insulated by multi-layer superinsulation and a vacuum 
space. The thermal vacuum jacket is supported by a 
PMMA acrylic vacuum chamber, 12 inches in diame- 
ter, as shown in Fig. [1] The inner face of the acrylic 
light guide rests against the window to the liquid argon 
cell. The light guides are 0-ring sealed to PMMA flanges 
on the insulating vacuum chamber, and the outer light 
guide face is at laboratory atmosphere and room temper- 
ature. The PMTs are coupled to the light guides with 
Bicron BC-630 optical gel and operated at room tem- 
perature. The light guides serve multiple purposes: to 
transport visible light from the liquid argon (at 87 K) 
to the room-temperature PMTs, as the vacuum seal for 
the insulating chamber, as a thermal barrier between the 
PMTs and the liquid argon, and to moderate neutrons 
from the PMTs. The measured heat load in this config- 
uration is approximately 7 watts per light guide. 

The target cell is filled with argon that has been pu- 
rified in the gas phase and then liquefied. The process 
system consists of a source of ultra-high purity (grade 
5) argon gas, a heated Zirconium metal SAES getter 
and a cooling coil that is immersed in liquid nitrogen. 
The argon gas at the output of the purification system 
contains less than 0.1 parts per billion for most impu- 
rities. No additional purification of the argon was per- 
formed after the detector was filled. The liquid nitrogen 
is maintained at a constant pressure of approximately 
33.5 psia to regulate the temperature of the cryostat. 
Argon is continuously recirculated through the cooling 
coil to remove the heat introduced through the acrylic 
light guides. The results presented here are for 67 days 
of operation, over which the light yield and detector re- 
sponse were stable (see Fig. [11]) . 

The target chamber is surrounded by a neutron and 



7-ray shield consisting of a minimum of two layers of 
ultra-pure water in 12-inch cubical polyethylene contain- 
ers, for a total water shielding thickness of 60 cm. A 
Geant4 0, § Monte-Carlo simulation showed that the 
water shield is sufficient to reduce nuclear recoil events in 
the liquid argon from external-source (a,n) neutrons to 
less than one per year. The data presented here are from 
operation of the detector at Queen's University without 
any overburden for cosmic-ray shielding. Further data 
are being collected underground at SNOLAB at a depth 
of 6000 meters water equivalent. These data will be re- 
ported later. At that depth the cosmic-ray muon flux 
passing through the detector is expected to be less than 
one per year, reduced from approximately 10 per second 
at Queen's. 
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FIG. 1; DEAP-1 detector configuration. Liquid argon is con- 
tained in a stainless steel target chamber, which includes an 
inner acrylic cylinder and a difi'use reflector. Scintillation 
light from the liquid argon is wavelength-shifted by TPB on 
the inner detector surface and is transmitted to PMTs op- 
erating at room temperature through PMMA acrylic light- 
guides. The PMTs (5-inch ETL 9390B) are held in place by 
spring-loaded polyethylene supports, which maintain contact 
between the PMTs and light guides. 



III. DATA AND ANALYSIS 
A. Electronics and Trigger Configurations 

Signals from the PMTs are sent to Phillips 778 ampli- 
fiers. One output of each amplifier is sent to a LeCroy 
digitizing oscilloscope that samples the waveforms at 
1 GHz for 10 fis, many times the 1.6 /is lifetime of 
the longest component of scintillation light. The sec- 
ond output of the amplifier is passed to a linear fan- 
out with one of the two outputs being used for discrim- 
inator triggering and the other for more sophisticated 
triggers. The "DEAP-1" trigger is a signal above dis- 
criminator threshold in both PMTs, corresponding to 
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approximately 1 keVcc, well below the region of interest 
for this study. The remaining output from the linear 
fan-out for each PMT is passed to a summing ampli- 
fier and then to a shaping amplifier with an integration 
time of 6 /is. This summed signal is passed to a single- 
channel analyzer (SCA) and used as a veto for events 
far larger in energy than the region of interest. Cut- 
ting these high-energy events allows us to increase the 
effective data-taking rate by a factor of approximately 
3. The DEAP-1 trigger with this veto imposed is called 
the "DEAP-1 SCA trigger". 

The response of DEAP-1 to electromagnetic interac- 
tions in the liquid argon is calibrated by a 10 /zCi encap- 
sulated ^^Na source placed outside of the argon chamber. 
The source emits a positron and a 1274 keV 7-ray. The 
positron annihilates in the source container producing 
back-to-back 511 keV 7-rays. One of the 511 keV 7- 
rays scatters in the argon, and we use two additional 
detectors to "tag" the backwards going 511 keV 7 ray 
and the additional 1274 keV 7 ray. The source is placed 
on the axis of an annular Nal detector with dimensions 
3.75" X 11.5" X 12" (inner diameter x outer diameter x 
length) to detect the 1274 keV 7-ray. The annulus is di- 
vided into 4 sections along its axis with a PMT on each 
section. The signals from the 4 sections of the annulus 
are independently amplified (Phillips 778) and discrim- 
inated (Phillips 705) with the threshold set just below 
the 1274 keV peak. The outputs of the discriminators 
are passed to a logical OR (Phillips 775 with coincidence 
level 1) to generate the "annulus trigger". The four an- 
nulus sections also have a higher-energy discriminator 
which is used as a veto on the annulus trigger to reduce 
backgrounds from cosmic rays. 

A cylindrical Nal crystal of dimension 3.25" x 3" (di- 
ameter X length) coupled to a PMT is placed at one 
end of the annulus. This PMT signal is amplified and 
passed to a single-channel amplifier (Ortec 420) which is 
centered on the 511 keV peak to detect the "backwards" 
going gamma ray. The "back PMT" trigger is thus gen- 
erated. The "global tag" is generated when an annulus 
trigger and a back PMT trigger are coincident. The ge- 
ometry of this triple-coincidence tag is shown in Fig. [2l 
The location of the source and back PMT are optimized 
to maximize the coincidence rate. 
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Recording with the oscilloscope is triggered by either 
the DEAP-1 trigger, the DEAP-1 SCA trigger or one 
of these two triggers in coincidence with the global tag. 
The timing of the triggers is set such that the leading 
edge of the trigger is defined by the back PMT so that 
approximately 1 /is of the waveform is recorded before 
the leading edge of the event. Pulse shape discrimination 
(PSD) calibration is carried out with the DEAP-1 SCA 
trigger in coincidence with the global tag. 
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FIG. 2: Calibration source geometry. The decay of ^^Na is 
tagged by 7's in the back Nal PMT, the annulus detector, 
and in DEAP-1, which allows a very low non-7 background 
for calibration of the pulse-shape discrimination. 



B. Data Flow 

The signal from each PMT is recorded on the os- 
cilloscope at a "high-gain" setting of 50 mV/division. 
When higher-energy signals are of interest the signals 
are recorded at both the high-gain setting and a low- 
gain setting of 500 mV/division on a second channel. 

A Linux data-acquisition computer runs an event 
builder that reads the files in real time and converts 
them to a ROOT-based ^ format. The root-based 
data files are automatically transferred to a disk farm 
on the High-Performance Computing Virtual Labora- 
tory (HPCVL) [10] for offline analysis. For the triple- 
coincidence calibration, the rate of scattering events in 
the argon is approximately 10 kHz, the rate of back-PMT 
and annulus coincidences is approximately 100 Hz, and 
the triple-coincidence rate is approximately 20 Hz. 



C. Analysis 

Figure [3] shows PMT traces from a sample 7-ray event 
in the region of interest from PSD calibration. A lin- 
ear baseline correction is applied waveform- by- wave form: 
the baseline is found in the first 0.5 /is of the waveform 
trace. That baseline is used to search for signals above 
threshold in the final 4 /is of the pulse. Areas around 
these signals are removed and a baseline is calculated at 
the end of the pulse. A linear function is found from the 
two baseline regions and subtracted from the total trace. 



The leading edge for each signal is found using a volt- 
age threshold of 5 mV, which corresponds to approxi- 
mately one tenth of a photoelectron. The prompt region 
of the pulse is defined from 50 ns before the edge to 150 



3 



> -50h 



Prompt region shaded 



PMT A 



T) 1 2 3 4 5 6 7 
Time (|J,s) 



9 10 



cussed later) are passed through a cut requiring \dt\ < 
20 ns. The average leading edge time of the two PMTs 
tedgo is used in PSD data to define a time window in soft- 
ware for accepting events in DEAP-1 with respect to the 
time of the back-PMT signal. For PSD data, a cut is im- 
posed that accepts only events whose leading edge time 
is within 30 ns of the mean time as found by a Gaussian 
fit to the distribution. This cut reduces random pileup 
of calibration 7 rays with high-Fpiomptbackgrounds, de- 
scribed later. Figure [4] shows the timing distributions 
for a typical PSD run. 
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FIG. 3: A sample event from a Na PSD calibration run. 
Shown are the recorded traces for each of the two PMTs, 
labeled A and B. The shaded area shows the prompt region 
which begins 50 ns before the leading edge of the pulse and 
ends 150 ns after the leading edge. 



ns after the edge. The late region is defined from 150 
ns to the end of the waveform, ~ 9 /is after the leading 
edge. The total charge for each PMT for both prompt 
and late regions is found and converted to units of photo- 
electrons using the single photoelectron charge calibra- 
tion described below. Defining Q*'f„,pt ^^^^ as the number 
of photoelectrons in PMT A or B that are prompt or late, 
and Total PE as the total number of photoelectrons, we 
write the fraction of prompt hght as 



Total PE 



(1) 



The relative signals in the two PMTs are used to re- 
construct the position of the event, Zgt, along the cylin- 
drical axis of the detector, 



35.2 cm X 



Total PE' 



(2) 



where Q* "" = Q^;^^^^ + Q*f and 35.2 cm is the distance 
from the center of the cylinder to the front face of the 
PMT. 

The time between leading edges St for the two PMTs 
is found. All data (PSD and neutron-source data, dis- 
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FIG. 4: The timing distributions for a typical run. The ver- 
tical lines indicate the cut values. For the top plot the cuts 
are always imposed at ± 20 ns. For the bottom plot the 
mean value is determined by a Gaussian fit and the cuts are 
imposed at the mean ± 30 ns. 



The mean charge from single photoelectrons (spe) is 
derived from these data. Using only pulses with a low 
number of photoelectrons (below 200), we scan the final 
4 /is of the pulses (high-gain channel only) looking for 
signals crossing a threshold of 5 mV. These signals are 
integrated from 5 ns before the threshold crossing to 20 
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ns after the threshold crossing (a total of 26 samples). 
On a run-by-run basis these spe charges are collected 
in a histogram. The mean charge is determined by the 
iterative procedure outlined in Ref. [ll| that finds the 
mean between 0.4 and 2.5 times the previous value until 
convergence is achieved. Figure [5] shows histograms for a 
typical run. The mean charges for the two PMTs are 1.73 
and 1.82 pC, with a variation throughout the running 
period of less than 4%. 
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FIG. 5: Single photoelectron spectra of the DEAP-1 PMTs, 
used to calibrate the absolute light yield and number of de- 
tected photons. 



The energy response is calibrated throughout 
the data taking with the full-energy peak from 
^^Na. The light yield is approximately 2.8 ± 0.1 
photoelectrons/keVee (see Fig. [TT]) where keVoc is the 
electron-equivalent energy, calibrated by the ^^Na peak. 
The light yield was stable throughout the 67-day running 
period. 



parameters that can be used to describe the Fprompt dis- 
tribution as a function of energy. As shown in SectionlVl 
we expect to use pulse-shape discrimination for signals of 
approximately 120 photoelectrons and above. We there- 
fore pre-determined the region of interest for demonstra- 
tion of PSD to be 120-240 photoelectrons, correspond- 
ing to approximately 43-86 keVoc in DEAP-1. Figure [5] 
shows the Fprompt distribution for 7-rays and nuclear re- 
coils, for events between 120 and 240 photoelectrons. We 
find that none of the 16.7 million 7-ray events leak into 
the nuclear recoil region, and infer that the PSD in liquid 
argon is less than 6 x 10~* from these data. Table IIVI 
shows the mean Fprompt values and the inferred sin- 
glet / triplet ratios for electronic and nuclear recoil events 
versus energy in liquid argon. 

We have also repeated this analysis with a lower photo- 
electron threshold and higher Fprompt region. For a neu- 
tron detection efficiency of 50%, there is one contamina- 
tion event observed between 70-240 photoelectrons (ap- 
proximately 25-86 keVoo), for an observed pulse-shape 
discrimination of 4.7 x 10^^ in that range. 
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FIG. 6: Fprompt versus energy distribution for the triple- 
coincidence 7-ray events. The region between 120-240 pho- 
toelectrons for Fprompt > 0.7 contains no events. 



IV. RESULTS 



Data were collected with the DEAP-1 detector be- 
tween August 20th and October 16th, 2007. A total 
of 63,072,900 triple-coincidence ^^Na calibration events 
were recorded for the PSD analysis. Figure [5] shows the 
distribution of Fprompt versus number of photoelectrons 
for these events; for this data set a hardware threshold 
imposed with the SCA restricts the high-photoelectron 
events to approximately 300 photoelectrons. Analysis of 
neutron calibration data shows that this does not remove 
high-Fprompt (0.7 < Fprompt< 1) uuclear recoil events. 
Figure [7] shows the Fprompt versus photoelectron dis- 
tribution for neutrons and 7-rays from an Americium- 
Beryllium (Am-Be) calibration source. From these data 
we have evaluated the PSD in liquid argon and measured 



Several tests of detector stability were performed. Im- 
purities in the argon can decrease the light yield by ab- 
sorbing UV photons or by quenching argon dimers before 
they decay, thereby also decreasing the observed lifetime 
of the triplet state in liquid argon |1^ . [l3l | . Therefore, we 
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FIG. 7: Fprompt versus energy distribution for neutrons and 
7 rays from an Am-Be calibration source. The upper band 
is from neutron-induced nuclear recoils in argon, the lower- 
band is from background 7-ray interactions. 



FIG. 9: Comparison of Zflt distribution for 7-rays from the 
PSD data, and for high-Fprompt backgrounds during the run 
(labeled Surface backgrounds). Also shown, for reference, is 
the distribution of high-Fprompt background events with the 
detector operating underground at SNOLAB. 
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FIG. 8: Fprompt distribution for 16.7 million tagged 7-ray 
events from the ^'^Na calibration, and nuclear recoil events 
from the Am-Be calibration, between 120 and 240 photoelec- 
trons (approximately 43-86 keVee). No 7-ray events are seen 
in the nuclear recoil region. 

measured the triplet lifetime in DEAP-1 over the course 
of the run to check that impurities did not build up in 
the detector over time. 

We use ^^Na calibration data to measure the triplet 
lifetime. For each calibration run, we find all events that 
pass the data cleaning cuts and contain over 200 photo- 
electrons. The raw traces for these events are aligned ac- 
cording to the measured trigger positions and summed. 
We then fit the following model to the average trace be- 
tween 500 and 3000 ns from the trigger: 

/(O = Aexp(-i/T3) + B, (3) 

where A is a normalization factor, ra is the triplet life- 
time and i? is a constant baseline term. 

As a consistency check, we measured T3 for photo- 
electron bins of size 200 between 200 and 1600 photo- 



FIG. 10; High-Fprompt background event rate versus time. 
The average background rate is 4.6 ± 0.2 mHz. 



electrons and did not observe any systematic effect from 
the signal size. There are systematic errors associated 
with both the fit window and the linear baseline correc- 
tion discussed in Section [III CI We estimated the size of 
the error associated with the fit window to be 40 ns by 
changing the start and end times of the fit by 500 ns. 
We performed the fit for both corrected and uncorrected 
traces and estimated the size of the error associated with 
the basehne to be 50 ns. We added the two estimated 
systematic errors to determine a combined systematic 
error of 60 ns. 

The measured lifetimes over the course of the run 
for traces without the baseline correction are shown in 
Fig-ttH in which the error bars shown are statistical only. 
We observe no significant increase in the impurity level 
throughout the run, and we measure the long time con- 
stant to be 1.46 lb 0.06 (sys) /iS, consistent with previous 
measurements [l^ [l^ • Further analysis of systematic 
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FIG. 11: Light yield stability during the run. The average 
light yield is 2.8 photoelectrons/keVcc- Shown are statistical 
uncertainties only. The scatter of the data points is within 
the systematic uncertainty of the measurement. 



uncertainties related to the fit window and baseline cor- 
rection as well as hardware improvements and investiga- 
tion of the effect on triplet lifetime of possible impurities 
in the argon are being investigated for the SNOLAB run. 
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FIG. 12: Measured lifetime of the triplet component in liquid 
argon during the run. Uncertainties shown are statistical 
only. The scatter in the data points is consistent with the 
systematic uncertainty of the measurement. No significant 
change in the lifetime was seen. 



TABLE 1: Measured values of Fprompt and singlet /triplet 
intensity ratios (/1//3) versus energy (T) for neutron-induced 
nuclear recoils and for electrons in liquid argon. 



electron recoils 



T 

(keVec) 


Fprompt 


I1/I3 


T 
(keV.) 


Fprompt 


I1/I3 


15 


0.38 


0.46 


50 


0.66 


1.71 


20 


0.36 


0.42 


66 


0.74 


2.50 


25 


0.35 


0.39 


83 


0.76 


2.80 


30 


0.34 


0.37 


100 


0.78 


3.17 


35 


0.33 


0.36 


116 


0.80 


3.49 


40 


0.33 


0.35 


133 


0.80 


3.66 


45 


0.32 


0.35 


150 


0.81 


3.79 


50 


0.32 


0.34 


166 


0.82 


4.05 


75 


0.32 


0.33 


250 


0.84 


4.57 


100 


0.31 


0.32 


333 


0.85 


4.88 


125 


0.31 


0.32 


416 


0.85 


5.01 


150 


0.31 


0.31 


500 


0.86 


5.37 


200 


0.29 


0.28 


666 


0.86 


5.56 



nuclear recoils 



PSD in argon. The average background rate in the re- 
gion of interest (120-240 photoelectrons) measured for 
the data presented in this paper is 4.6 ± 0.2 mHz, con- 
stant throughout the run as shown in Fig. [101 

With the triple-coincidence tag, the measurement 
on surface is limited by random coincidences between 
the global tag (back-PMT and annulus) and high- 
Fprompt background events in DEAP-1. The rate of 
high-Fprompt pile-up during the 7-ray calibration run is 



R 



bkg 



RtagRn/^t, 



(4) 



where i?tag is the global tag rate, i?„ is the rate of high- 
Fprompt background events measured in DEAP-1 and 
At is the width of the time window (applied with the 
software cut). 

The achievable PSD can be found by dividing i?bkg 
by the rate of data acquisition in the energy region of 
interest (ROI), i?Roi: 



D 



bkg 



-Rbkg _ RtagRn^t 



Rroi 



Rroi 



(5) 



High-Fprompt (nuclear recoil) backgrounds are evalu- 
ated using the DEAP-1 trigger without any calibration 
sources present. Figure [5] shows the Z^t distribution 
of high-Fprompt background events, compared to that 
of ^^Na calibration 7-ray events. Also shown for ref- 
erence is the Zfit distribution of high-Fprompt back- 
ground events with the DEAP-1 detector operating un- 
derground at SNOLAB. The reduced backgrounds un- 
derground will allow a more sensitive measurement of 



Dbkg is thus the discrimination level where we would 
expect to find one background event. Table |TT] shows 
the relevant parameters and the result from calculating 
Dbkg, and this background is below the demonstrated 
discrimination level of 6 x 10^*. 

It should be noted that the high-Fprompt background 
rate has been reduced by a factor of approximately 10 at 
SNOLAB, and further background reduction is expected, 
allowing for a more sensitive measurement of PSD in 
argon. 
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TABLE II: Background from random coincidence between 
the global tag and high-Fprompt background events in the 
argon. -Rtag is the coincidence rate between the annulus and 
back-PMT with the ^^Na source in place, At is the coinci- 
dence time window imposed by the analysis, Rn is the rate of 
high-Fprompt background events in the liquid argon, -Rroi is 
the rate of triple-coincidence events and Dbkg is the PSD level 
where one radon-coincidence background event is expected. 



we have 



Variable 


Values 


iJtag (Hz) 


1000 


At (ns) 


60 


R„ (mHz) 


4.6 


Rnoi (Hz) 


18 


-Dbkg 


1.53x10"* 



V. ANALYTIC MODEL FOR 
DISCRIMINATION POWER 

We developed a model for the Fpj-ompt distribution by 
assuming that the singlet and triplet contributions are 
both Gaussian and uncorrelated. The distribution of 
Fprompt values Fp is then a ratio of two Gaussian distri- 
butions, which to a good approximation can be written 
as [11 



1 



'2tt {erf 



(6) 



The distribution is thus a function of four parameters: 
the means /ip, hl and widths crp, gl of both the singlet 
(prompt light) and triplet (late light) distributions. Each 
of the four parameters depends on the event energy, and 
the ideal distribution in Eq. ([6]) has to be convolved with 
the energy response of the detector to adequately model 
the observed data. 

By definition, the parameters fip,^L are functions of 
the most likely Fpiompt value Fp-. 



MP 



F N 



and 



/ii = (1 - Fp)Npe 



(7) 



(8) 



so that their energy dependence is known if that of Fp 
is known. 

The energy dependence of the distribution widths 
<^Pt<^l is modeled by splitting a up into the following 
components: uncertainty in the single photoclectron cal- 
ibration CTspe, uncertainty in the counting window due to 
uncertainty in the trigger time (Twin , electronic noise CTcIcc 
and statistical uncertainty (Tstat , so that 



spc 



clcc 



(9) 



Writing the statistical noise as the square root of the 
number of photoelectrons in the prompt or late window 



2 

"stat.P 
<^stat,L 



F N 
-'- p^ ^ pe 



(1 - Fp)Npe. 

The single photoclectron noise contribution is then 



(10) 
(11) 



c^spcP = ia^p^/^j.^p^)jFpNpe 



^spc,L 



«c/Mspc)\/(l-^p)^P- 



(12) 

(13) 



where /i^^ and cr^j, are the mean and width of the single 
photoclectron peak. 

The uncertainty in the counting window of approxi- 
mately 30 ns leads to an uncertainty in the late photo- 
electron count of 2.5%. Electronic noise of 0.58 photo- 
electrons in the prompt window, and 3.7 photoelectrons 
in the late window are measured by triggering the de- 
tector with a pulser and analyzing the PMT traces. The 
fraction o-spc/l^spc measured to be 0.276. The final ex- 
pressions for the dependence of cp and ctl on the number 
of photoelectrons, and thus on the energy of the event, 
are then 



al ^FpNpe + {0.276^J FpNpe 
+ (0.025(1 - Fp)Npe)^ 



(0.58 pe)^ 



= (1 - Fp)Np, + (0.27601 -i?p)iVpe)' 
+ (3.7 pe)2 + (0.025(1 - Fp)Npef 



(14) 



(15) 



The behaviour of the prompt and late noise versus the 
number of photoelectrons is illustrated in Fig. [13] and 
Fig.m 

Finally, since the Fpiompt distribution depends on en- 
ergy and number of photoelectrons, the distribution is 
convolved with the energy spectrum of interest N{E): 



P*{Fp,Np,) = Y,N{n)P{fp{n),, 



n=0 



cr(n) v27r 



(16) 



where the energy E has been converted to number of 
photoelectrons n and a(E) is the energy resolution func- 
tion as determined from the detector's energy calibra- 
tion. This model agrees with the data for energies above 
60 photoelectrons. A comparison between the model and 
the data is shown in Fig.[Tn]for the energy region 60-120 
photoelectrons, and in Fig. [16] for the energy region 120- 
240 photoelectrons. The dark gray line is from Eq. 
and the shaded area represents the systematic uncer- 
tainty from uncertainties in energy calibration and noise 
parameters. 
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FIG. 13: The noise level ap in the prompt time window is 
shown as a function of the total number of photoelectrons 
in the pulse. The broken lines show how the four sources 
of noise and uncertainty that are taken into account in our 
model contribute to the total noise level. 
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FIG. 14: The noise level ctz, in the late time window is shown 
as a function of the total number of photoelectrons in the 
pulse. The broken lines show how the four sources of noise 
and uncertainty that are taken into account in our model 
contribute to the total noise level. 



FIG. 15: Comparison of the data and the analytic model 
in the region 60-120 photoelectrons (approximately 21- 
43 keVcc). The gray shaded area represents the uncertainty 
in the noise parameters and in the energy calibration. 



c 

O 
U 




0.2 



prompt 



0.4 



0.6 



0.8 



projection at 120 to 240 pe 



FIG. 16: Comparison of the data and the analytic model 
in the region 120-240 photoelectrons (approximately 43- 
86 keVeo). The gray shaded area represents the uncertainty 
in the noise parameters and in the energy calibration. 



parameters are negligible, the acceptance for nuclear re- 
coil events, and the limit expected from surface back- 
ground contamination. 



To represent the discrimination power graphically, we 
define the variable Picak as 

PlcukiF') = /o Fprompt(a:)(ix ^^^^ 

for a given Fprompt distribution. The value of Picak for a 
given Fpi-ompt value thus represents the pulse-shape dis- 
crimination. Figure [TT] shows the Pioak distribution de- 
rived from the ^^Na calibration data, along with the pre- 
dictions from the analytic PSD model. Also shown are 
the analytic model assuming that the systematic noise 



VI. DARK MATTER SENSITIVITY 

Using Eq. ^6]) we next evaluate the expected discrim- 
ination power in liquid argon versus energy threshold as- 
suming a large target mass detector. We have designed 
DEAP-3600, a large spherical detector consisting of ap- 
proximately 250 PMTs surrounding a spherical target 
with a mass of 3600 kg of liquid argon. Gcant4 simula- 
tions benchmarked against the light yield in the DEAP-1 
detector predict that in this geometry a light yield of ap- 
proximately 6 photoelectrons/keVcc could be achieved. 
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FIG. 17: Picak distribution from Na calibration data from 
DEAP-1, and analytic models with and -without additional 
noise parameters for 120-240 photoelectrons. The lo-wer curve 
shows the expected backgrounds in the measurements from 
high-Fprompt events. 



FIG. 18: Pioak distribution expected for 20-40 keVcc with 
6 photoelectrons/keVco yield. Shown are the distribu- 
tions with and without systematic noise terms, and the 
Fprompt values for 50% and 90% nuclear recoil acceptance. 



Assuming the detector is constructed of clean materials 
and appropriately shielded so that genuine nuclear recoil 
backgrounds have been mitigated, the dominant detec- 
tor background in argon -will be from (3 decays of "^^Ar. 
Argon that is condensed from the atmosphere is kno-wn 
to contain cosmogenically-produced '^^Ar, with a rate of 
decays of approximately 1 Bq per kg jl^ . Recent -work 
has sho-wn that there is a possibility of obtaining argon 
that has been sequestered underground that is depleted 
in '^^Ar by a factor of 20 [l3|- An optimization analysis 
has sho-wn that we could achieve an analysis threshold 
of approximately 20 keVee for a dark matter particle 
search -with natural argon; from Eq. p6)l the expected 
discrimination is 10"^*^ at that threshold, -which is suffi- 
cient to reduce the backgrounds from atmospheric '^^Ar 
decay to less than 0.2 events in a three-year run with 
1000 kg of liquid argon in the energy window from 20- 
30 keVco- In that analysis we have assumed that for 
the large detector, each of the approximately 250 PMTs 
would detect on average less than a single photoelectron, 
and thus that the systematic noise parameters can be 
reduced so that they are negligible compared to the sta- 
tistical uncertainties in Eq. The Picak distribution 
given by Eq. ([T7]) is shown in Fig. [TH] Fig. [TO] shows the 
90% confidence level sensitivity versus WIMP mass for 
both cases, compared to the current experimental lim- 
its of the CDMS [H] and XENON-10 {m] experiments. 
We have used standard assumptions for the galactic halo 
outlined in reference [201 and assumed a nuclear recoil 



scintillation efficiency of 30% compared to electrons [21 1 
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FIG. 19: Dark matter sensitivity of liquid argon. Shown 
are the current experimental limits from the CDMS and 
XENON-10 collaborations, and the expected sensitivity for 
1000 kg of liquid argon with a 20 keVee threshold, and with 
a 19 keVec threshold for argon that has been depleted in 
^''Ar by a factor of 20. 



VII. SUMMARY 

We have shown that pulse-shape discrimination in liq- 
uid argon can be used to separate 7-ray and nuclear 
recoil events in the energy region 43-86 keVco with a 
leakage of less than 6 x 10^^ for full recoil detection 
efficiency, and in the energy region of approximately 25- 
86 keVoo with a leakage of less than 4.7 x 10~* for a 
50% recoil detection efficiency. An analytic model which 
assumes simple Gaussian distributions for the singlet 
and triplet components of the scintillation signal and 
accounts for statistical and systematic sources of noise 
in these distributions agrees with these data for events 
above 21 keVoc. This model projects that available dis- 
crimination power in argon is 10~^'^ at 20 keVcc, for a de- 
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tector with 6 photoelectrons/keVee light yield, allowing 
for sufficient background rejection of '^^Ar in a 1000-kg 
liquid argon dark matter experiment. 
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